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ABSTRACT 

I study some aspects of common envelope evolution, where a compact star enters the envelope 
of a giant star. I show that in some binary systems under a narrow range of parameters, a substantial 
fraction of the giant stellar envelope is lost before the onset of the common envelope. The reduced 
envelope mass at the onset of the common envelope implies that the binary system emerges from the 
common envelope with a relatively large orbital separation. I therefore caution against a simple 
treatment, which omits this process, in the study of systems that eveloved through a common 
envelope phase and ended with a relatively large orital separation, e.g., PG1115+166 . A fraction 
of the envelope that is lost while the companion is still outside the giant stellar envelope is accreted 
by the companion. The companion may form an accretion disk and blow two jets. I propose this 
scenario for the formation of the bipolar planetary nebula NGC 2346, which has a binary nucleus 
with an orbital period longer than that of any other known binary system in planetary nebulae. 
Because this scenario, of a late common envelope phase, occurs for a narrow range of parameters, 
I expect it to be applicable to a small, but non-negligible, number of systems. 

Subject headings: stars: binaries — stars: evolution — stars: AGB and post-AGB — stars: 
mass loss 
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1. INTRODUCTION 



Many evolved close binary systems acquire their close orbit through common envelope evolution 
(Paczynski 1971; Iben & Livio 1993; Armitage & Livio 2000; Taam & Sandquist 2000). In the 
common envelope phase a compact companion enters the envelope of a more extended star, e.g., 
an asymptotic giant branch (AGB) star, and because of tidal interaction and friction, the orbit 
shrinks (see review by Iben & Livio 1993). In principle, orbital energy and angular momentum are 
transferred to the giant stellar envelope, where part of the orbital energy goes to expel the envelope. 
A commonly used parameter is the ratio of the orbital energy that is released during the common 
envelope phase, AE^t,, to the binding energy of the ejected envelope, AE^nd'- «ce = A.Ebmd/ A-E rb 
(e.g., Livio & Soker 1988). Note that different definitions for the binding energy exist (e.g., O'Brien, 
Bond, & Sion 2001). Since the orbital energy that is released depends mostly on the final orbital 
separation, the value of «ce can in principle be calculated for systems whose final orbital separation 
is known. This value can then be used to derive the final separation in the study of the common 
envelope phase in binary systems evolution. To calculate the value of «ce in a specific system, the 
giant structure at the onset of the common envelope should be known, in particular the envelope 
mass and radius. Using different values for the envelope mass and radius yields a very different 
value for a C E (O'Brien et al. 2001; Maxted et al. 2002). 

Maxted et al. (2002) assume that negligible mass has been lost prior to the onset of the 
common envelope in PG1115+166, and find «ce > 1> so they propose, that internal energy stored 
in the envelope, e.g., ionized gas, contributes to envelope ejection. I disagree with this proposal 
for the following reasons. First, Harpaz (1998) convincingly argues that the energy released in the 
recombination of the envelope does not contribute much to the ejection of the envelope. Second, 
on a more general ground, energy is not a problem in the mass loss process of AGB stars. In 
most AGB stars, the momentum carried by the wind, M w v w , is not much above the momentum 
supplied by the radiation L/c, where M w is the mass loss rate into the wind, v w the wind speed, 
c the speed of light, and L the stellar luminosity. In most systems M e v w <J 3L/c (Knapp 1986), 
implying that the ratio of stellar luminosity to the rate of kinetic energy carried by the wind is 
L/E w ^ 10 4 , assuming v w ~ 20 km s -1 . Thus, the energy is not a problem as regards the ejection 
of an AGB stellar envelope, and it is not clear how relevant the ckec parameter is by itself. This 
can be demonstrated for PG1115+166, for which the final masses are O.43M and O.52M , and the 
final orbital separation is a/ = 45-R© (Maxted et al. 2002). The orbital energy that was released 
during the spiraling-in process was A-E^b — 10 46 erg. An AGB star with a luminosity of 5000-L Q 
releases this energy in only < 20 years. 

In any case, another way to reduce the calculated value of acE is to reduce the binding energy 
of the envelope by considering a decrease of its mass before the onset of the common envelope. This 
is the goal of the present paper. I show that for a small initial parameters space, the companion 
brings the giant envelope to corotate with the orbital motion, hence the spiraling-in process stops 
before the companion reaches the envelope, such that only after the giant star loses a substantial 
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fraction of its envelope does the compact companion enter the envelope. In §2 I discuss the initial 
stage of reaching synchronization, and in §3 I discuss the onset of the common envelope. The 
general evolution of a close companion to a giant star was treated before, e.g., Alexander, Chau, 
& Henriksen (1976), who include accretion by the companion, drag, tidal force, and eccentricity. 
They note that the effect of a wind blown by a rotating giant is to reduce the orbital separations, 
but do not present a detailed solution for the orbital separation. In the present paper I neglect 
some effects studied in the past. This allows a simple solution for the orbital separation, which 
clearly demonstrates the points raised here. I discuss this solution and its application to specific 
systems in §3. My short summary is in §4. 



2. REACHING SYNCHRONIZATION 



I treat systems in which the secondary compact star, i.e., a main sequence star, a white dwarf, 
or a neutron star, is massive enough to bring the mass-losing giant star to synchronization. I neglect 
eccentricity, mass accretion and accretion drag by the companion, effects which were studied in the 
past (e.g., Alexander et al. 1976). I do consider tidal interaction, assuming synchronization is 
reached, and the effect of mass loss from a rotating giant. The system is taken to be stable to 
the Darwin instability. This stability criterion is I or b > 3/ e nv (e.g., Eggleton & Kiseleva-Eggleton 
2001), where I or b is the moment of inertia of the binary system, I cnv = r]M env R 2 is the moment of 
inertia of the giant, R g is the giant's radius, M cnv is the envelope mass, and r/ ~ 0.2 for AGB stars. 
After substituting for the orbital moment of inertia, the criterion for Darwin stability reads 

M » M 2 ..2 ,„.., „2 



where a is the orbital separation, and M g the giant stellar mass. Because of tidal interaction, orbital 
angular momentum is transferred to the giant's envelope, and the orbit starts to shrink when the 
giant's radius R g increases to R g ~ 0.1 — 0.2a. The exact ratio R g /a when synchronization and 
circularization occur strongly depend on the eccentricity (Hut 1982; Eggleton & Kiseleva-Eggleton 
2001), and weakly on some of the giant's properties, e.g., evolution time (see scaled equations 
in Soker 1996). After the system reaches synchronization, and the spiral-in process stops for a 
while, the giant radius is likely to further increase. I take an orbit with initial semi-major axis 
of Oj, while the radius when synchronization and circularization are achieved is ao < a^. To 
derive a simple expression, I take M<i <C M g at this evolutionary phase, before the giant loses 
most of its mass. The orbital angular momentum released in the spiraling-in from ai to ao is 
AJorb — M2{GM g ) l l 2 (a l J 2 — Oq 2 ). This equals the angular momentum gained by the envelope, 
which, assuming that initially the envelope has no angular momentum, is J cnv = I c11v ujo, where 
ojq = (GMg/afy 1 / 2 is the angular velocity of the envelope at the onset of synchronization. The 
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equality AJ orb = J cnv gives 



ao) 



1/2 



R„ 



0.2M C] 
M 2 



(2) 



Strong tidal interaction starts at ao ~ 5 — 7R g (Soker 1996). For M 2 = 0.1M cnv and ai = 7R g , 
synchronization is reached at ao ~ 6R g by equation (2). The system is still stable to the Darwin 
instability (eq. 1). Much lower mass secondary stars will spiral-in all the way to the envelope, 
because ao will be smaller and the orbit becomes unstable by equation (1). For the envelope to 
influence the orbital motion, its mass can't be much below that of the companion. Assuming that 
most of the initial giant mass M g is in the envelope, I find that the arguments of this paper arc 
applicable to secondary stars in the mass range 

0.1M 90 <S M 2 £ 0.5M 90 , (3) 

where M g o is the initial giant mass (when synchronization occurs). 



3. ENTERING A COMMON ENVELOPE PHASE 



Assuming again a synchronized circular orbit, the angular velocity of both the orbital motion 
and the giant's envelope is 

"G(M 2 + M 9 ^ 1/2 



LO 



(4) 



The wind blown by the giant carries angular momentum at a rate given by 



iwmd = (a 2 g + [3R 2 a)vM w , 



(5) 



where a dot stands for a time derivative, M u 



-M„ is the rate of mass blown in the wind, 



a g = aM2/(M g + M 2 ) is the distance of the mass-losing giant star from the center of mass, and (5 
is a factor which depends on the mass-loss geometry: (3 = 2/3 for a uniform mass loss from the 
giant's surface, while (3 = 1 for a mass loss from the equatorial plane. The first term inside the 
parenthesis of equation (5) is the angular momentum of the wind due to the motion around the 
center of mass, while the second term is due to the rotation of the envelope. The sources of the 
angular momentum carried by the wind are the orbital and giant envelope's angular momentum. 
The orbital angular momentum changes at a rate of 



orb 



M„ 



Jorb M g 



M g a 
2(M g + M 2 ) + 2a"' 



(6) 
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where J orD is the orbital angular momentum, and, as earlier, I assume a circular orbit. Neglecting 
the change in R g and in the core mass, M core = M g — M cnv , the rate of change of the giant's 
envelope angular momentum is 

where r\ is defined before equation (1). Conservation of angular momentum implies the equality 
•/orb + Jenv + Jw = 0. Using equations (5)-(7), the last equality gives the expression for the rate of 
change of orbital separation a as 



M g 



Mg + M 2 

(M g + M 2 ) 2 (R g \ 2 (M g + M 2 ) 2 (R g \ 2 M g + M 2 (R g \ 2 

2/J M„M 2 {a J +2V M„M 2 {a J + vMmv M q M 2 {a J 



(8) 



For a non-rotating giant and/or R g <C a, the familiar expression (a/a) = —M g /(M g + M 2 ) is 
recovered, and the orbit expands with mass loss. For close orbits, where both R g is not much 
smaller than a and the giant rotates fast, the second term in the square brackets on the right hand 
side dominates, and the orbit shrinks. The binary system may enter a common envelope phase, 
but only after some fraction of the envelope mass has been lost. 

To cast equation (8) in a simple form, I define 

X = l-3r ? M env » . (9) 



M g M 2 V a 

The stability against Darwin instability, equation (1), is x > 0. I also define 

Because the accretion by the companion is neglected, M g = M 2 fi. Canceling the time derivative, 
equation (8) can be written in the form 

x ^ = _^ + 4 ^fci)_ 4 ^ii±i_ 2 ^^I. („) 
a/! n + l L 1 A 4 M g 

For the effects studied here to be important, 2 ^ fi ^ 10 (eq. 3). Taking this range for the value 
of n, and the range for the the value of the wind geometrical factor (2/3) < (3 < 1, I approximate 
l)//x ~ 1. In light of the effects neglected here, e.g., accretion by the companion, I also make 
the approximations that x is constant, fi + 1/n ~ 1, and M env /M g ~ 1. These assumptions yield a 
simple form for equation (11): 

da 2 2a 2 „ . „ 9 . , 
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which has a simple solution 

( a\ 2 (a<,y ( uo + l \Vx 4-6,, , ( fn + l \ 2 '" 4-6g 

(^J "UJ Utt) -TT7 ( "° + 1) l— j + TT7 (f, + 1) ' (13> 

where the initial condition, after synchronization and circularization are reached via strong tidal 
interaction, is a(//o) = ao- 



4. DISCUSSION 



When <C ao the first term on the right hand side of equation (13) dominates, % = lj and 
the familiar expression for orbit expansion with mass loss is recovered. For a general case, when 
condition (1) holds and < % < 1, then the factor [(/xo + l)/(/i + l)] 2 ^ x that multiplies the first two 
terms on the right hand side of equation (13) becomes very large compared with the third term. 
In that case, there are three cases which are determined by the ratio of the two first terms on the 
right hand side 

r 

If r = 1, the first two terms cancel each other, and the orbit shrinks according to ao — 
R g (pL + l) 1 / 2 , where I assumed (4 — 677) / (2 + %) ~ 1. A common envelope will occur only if the 
giant radius increases a little during this process, and V becomes > 1. In the first case, occurring 
for T < 1, there will be no common envelope phase. In the second case, for V S> 1 (as I show below, 
practically for T ^ 1.3), the orbit decreases fast with mass loss, and a common envelope starts after 
only a small amount of the envelope mass has been lost. 

In the third case, the one this paper aims at, the onset of the common envelope occurs after 
a large fraction of the envelope mass has been lost. Consider the following case, which I choose to 
demonstrate the third case. I take ao/R g = 2.5, f? = 0.2, jjlq = M g /M2 = 6, and an average, over 
time, value of % = 0.5. The value of ^lq was chosen to match that of the progenitor of PG1115+166 
(Maxted et al. 2002). These values give T = 1.25. The orbit shrinks, and the common envelope 
phase starts, i.e., a = R g , when /j, ~ 4.3. If the giant's core mass is M corc ~ M2, these values 
of fio = 6 and \i = 4.3 mean that the envelope mass decreased to ~ 2/3 its initial mass. For 
higher values of T, the common envelope phase will start at earlier stages (depending on the other 
parameters as well), with a small decrease in envelope mass. I conclude that the scenario I am 
aiming at, i.e., that ;> 50% of the envelope is lost before the onset of the common envelope, is 
applicable for a narrow range of parameters, for which 1 < T ^ 1.3. Because x> an d to lesser 
degree R g , change during the evolution, the exact range for T should be determined by numerical 
evolutionary calculations. I propose that the progenitor of PG1115+166 had parameters similar to 
the values chosen above, and in consequence a large fraction of the envelope mass was lost before 
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the onset of the common envelope. The low envelope mass at the onset of the common envelope 
means a low binding energy. This, I suggest, accounts for the relatively large orbital separation 
found by Maxted et al. (2002). 

The evolutionary scenario studied here is relevant also to the formation of some bipolar plane- 
tary nebulae (PNs). Bipolar PNs are PNs with two lobes and an equatorial waist between them are 
similar in many respects to bipolar nebulae around symbiotic systems (Corradi & Schwarz 1995; 
Corradi 1995), and are thought to be formed from binary systems. In most cases the binary com- 
panion to the AGB star is outside the envelope (Soker & Rappaport 2000). From the 16 PNs known 
to have close-binary nuclei (Bond 2000), only NGC 2346 is a bipolar PNs. The orbital period of the 
central-binary system in NGC 2346 is P Q rb = 16 day, much longer than those of the other 12 PNs 
with determined orbital periods, all having P or b < 3 day. It is likely, therefore, that the companion 
to the AGB progenitor of NGC 2346 entered the envelope when the envelope mass was already 
small, so it did not spiral-in as close to the AGB core as in the other systems. I propose that during 
most of the evolution time the companion was outside, but close to, the envelope, where it accreted 
mass and blown collimated fast wind (CFW), or jets, which shaped the bipolar lobes (Soker & 
Rappaport 2000). Because of the narrow range required for the above scenario, 1 < T <; 1.3, I do 
not expect to find many bipolar PNs with close binary systems at their centers. 



5. SUMMARY 



The main purpose of this paper is to caution against a simple use of the «ce parameter when 
studying the final separation of binary systems emerging from a common envelope. acE is defined 
as the ratio of the binding energy of the envelope to the orbital energy released by the spiraling-in 
binary system. 

There is no need to invoke an extra energy source if a simple calculation shows that in a specific 
system «ce > 1- First, plenty of energy is supplied by the nuclear burning in the giant's core, ^ 10 4 
times the rate of kinetic energy carried by the wind. The gravitational energy released by the 
companion spiraling-in inside the giant envelope, though, has the advantage of being transferred 
in part directly into mechanical energy, by spinning-up the envelope, by tidal interaction, and by 
exciting pressure modes. Possible jets blown by the companion also have this advantage (Armitage 
& Livio 2000). (A small fraction of the nuclear energy is also transferred to mechanical energy via 
stellar pulsation, which becomes very strong on the upper AGB.) 

Second, as I showed in this paper, in some systems a substantial fraction of the envelope may 
be lost before the onset of the common envelope, thereby reducing the envelope's binding energy, 
hence the value of the calculated «ce- The evolution is as follows. As the giant radius increases, 
tidal interaction becomes stronger. When the ratio of the giant radius to the orbital separation 
becomes R g /a ^ 0.1 — 0.2, and the giant mass to companion mass ratio is = M g /M.2 <^ 10, the 
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system reaches synchronization and circularization. The spiraling-in process stops. As mass loss 
proceeds, angular momentum is lost from the system. The giant most likely continues to expand 
and the ratio a/R g decreases, so the effective value of ao/R g becomes smaller than that at the 
moment synchronization is achieved. The condition for this scenario to occur is 1 < V ^ 1.3, where 
T is given in equation (14). If V < 1 the system does not enter a common envelope phase, while 
if T ^ 1.3 the system enters a common envelope phase after only a small amount of envelope has 
been lost. The exact condition of T depends on the values of the different parameters, as well as 
on some effects not studied here, e.g., accretion by the companion, drag, and eccentricity. 

This scenario is applicable to systems which went through a common envelope and end with 
a relatively large orbital separation, e.g., PG1115+166 (Maxted et al. 2002). When outside the 
envelope, the companion can accrete from the AGB wind, form an accretion disk and blow jets. 
This will form a bipolar planetary nebula. Later the companion enters the envelope and spirals-in, 
ending at closer orbit to the core. I proposed this scenario to account for the bipolar planetary 
nebula NGC 2346, which has a binary nucleus with orbital period of 16 days (Bond 2000). Because 
this scenario occurs for a narrow range of parameters, I expect it to be applicable to a small, but 
non-negligible, number of systems. 
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